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ABSTRACT 

The electron mean free path in the intracluster medium (ICM) of galaxy clusters is much larger 
than the gyroradius; thus, heat is transported anisotropically along magnetic field lines. We show that 
the intracluster medium is unstable to the magnetothermal instability (MTI) using MHD simulations 
with anisotropic thermal conduction. As a result of the MTI, we find that the temperature profile 
of the ICM can be substantially modified on timescales of several billion years while the magnetic 
field is amplified by dynamo action to more than fifty times the original energy. We also show that 
the instability drives field lines to become preferentially radial leading to conduction that is a highly 
efficient fraction of the Spitzer conductivity. As such, we present the first self-consistent calculation 
of the effective thermal conductivity in the ICM. 

Subject headings: convection — galaxies: clusters: general — instabilities — MHD — plasmas — X-rays: 
galaxies: clusters 
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1. INTRODUCTION 

Clusters of galaxies are the most massive 
gravitationally-bound objects in the universe, yet a 
mere 3% of their mass is in the form of stars. The 
majority of their mass (84%) is in dark matter with the 
remaining 13% consisting of a hot, low density, mag- 
netized plasma called the intracluster medium (ICM). 
Due to their enormous masses, surveys of clusters 
provide key tests of cosmology. With the Chandra X-ray 
Observatory, we are capable of measuring the spatially 
resolved X-ray luminosity (I0 43 -I0 46 erg/s) emitted 
by the ICM to determine its density as a function of 
radius. In order to fully understand the fundamental 
astrophysics of the largest bound objects in the universe, 
we must develop an accurate quantitative picture of 
the heating and cooling mechanisms as well as the 
dynamical processes that govern the ICM. Thus, we are 
motivated by the use of clusters as a cosmological tool 
to study the physics of the ICM much more rigorously. 

The intracluster medium is heated by gravitational in- 
fall such that typical temperatures are 1-15 keV and 
typical densities are in t he range of 10 _3 -10 -2 cm -3 
(Peterson fc Fabianl [20061 ) . Since the ICM is fully ion- 
ized and radiative forces are negligible, we can treat ra- 
diation with optically thin cooling and neglect electri- 
cal resistivity. Estimates for the magnetic field strength 
in clusters range from roughly 1-10 /zG at the center 
and 0.1-1.0 ^G at a radius of 1 Mpc, values that cor- 
respond to a plasma beta, (3 = 8irP/B 2 200-2000 a 
dynamically weak magnetic field (jCarilli fc Tavlorf 2002) . 
In this dilute plasma, however, the mean free path for 
electron collisions can be twelve orders of magnitud e 
larger than the gyroradius (|Naravan fc Medvedev 2001). 
In this circumstance, the equations of MHD must be 
supplemented with anisotropic terms that include the 
near free-streaming motions of particles along magnetic 
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field lines ([Braginskiil [1965) . The parallel thermal con- 
ductivity of the electrons is larger than that of the 
ions by a factor proportional to (mi/m,,) 1 / 2 , whereas 
the parallel viscosity of the ions is larger t h an th at of 
the electrons by the same factor. [Balbusl ( |2000[ ) has 
predicted that gravitationally bound plasmas in this 
regime are susceptible to the magnetothermal insta- 
bility (MTI). This convective instability is driven by 
anisotropic thermal conduction along field lines with 
major d ynamical consequences ([Parrish fc Stond I2005L 
200730; iChandran fc~Denn is 2006). Thus, we are no 
longer justified in considering the weak magnetic field in 
the ICM as a purely passive component in a fluid gov- 
erned solely by hydrodynamics. 

There are two long-standing puzzles about the ICM. 
First, the radial temperature pr ofiles show a sligh t de- 
cline in temperature with radius (|Pratt et al.ll2007f ) that 
is smaller than th at expected from structure formation 
calculations, e.g. (|Loken et al J 12002). Second, the mag- 
netic field has been significantly amplified in the clus- 
ter from its primordial value of < 1 nG. The ICM is 
unstable to the magnetothermal instability on scales of 
tens of kpc and larger, as calculated from the current 
observed magnetic field strengths. In three-dimensional 
simulations, we have shown that the MTI is capable of 
generating convective motions and a magnetic dynamo 
while efficiently transporting heat. In addition, the MTI 
can make the temperature profile more isothermal as it 
exhausts its source of free energy. These similarities pro- 
vide a strong motivation for studying the ICM with this 
type of simulation. 

Hydrostatic equilibrium dictates that 



M{r) = - 



G/irrip 



T e (r) 



dn e (r) 
dr 



n e (r) 



dT e (r) 
dr 



(1) 



wher e \i is the mean molecular weight (|LaRoque et al.l 
2006). When an X-ray observation is made, the den- 
sity profile (first term of eqn. [Tj) can be constrained 
fairly well due to the property that L x oc n 2 . for free- 
free emission, the dominant cooling mechanism. Often, 
the gas density profile is fit with a n isothermal-/? model 
(|Cavaliere fc Fusco-Femianolll976l ). On the other hand, 
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the temperature profile (second term of eqn. Q]) can 
only be constrained by spectra, which are increasingly 
difficult to measure at larger redshifts. For high redshift 
clusters, the virial masses are often constrained simply by 
assuming isothermality, in contradiction to local observa- 
tions. Why do we care about determing mass with such 
accuracy? A flux-limited X-ray survey c an be used to 
constr ain f^o, the matter parameter, e.g. iBorgani et al.l 
(1999). Locally, however, there is a degeneracy between 
the parameters fio, the matter fraction of the universe, 
and as, the variance of the power spectrum on 8/i _1 Mpc 
scales. Thus, one must examine the evolution of the clus- 
ter mass function with redshift to break this degeneracy. 

In iJ21 we review the physics of the MTI then procede 
to construct a model of the ICM for numerical simu- 
lation in fj3j In ||4] we outline our numerical methods. 
We then show in S}5] that the MTI is capable of rearrang- 
ing the temperature profiles of clusters on cosmologically 
important timescales, i.e. 0.5-10 billion years. We also 
show that a magnetic dynamo operates as a result of the 
turbulent convection driven by the MTI. A comparison 
to runs with purely isotropic conduction at a specified 
fraction of the Spitzer conductivity show that much less 
turbulence and no dynamo result. Finally, we point out 
that the magnetic field is driven to a highly radial geom- 
etry and speculate on possible observational implications 
before concluding. 

2. PHYSICS OF THE MTI 

The physics of the magnetot herm al instabi li ty ar e 
clearly outlined in iBalbusI pOOOh and IQuataertl (|2008f ). 
We will only briefly review them here. The magnetother- 
mal instability is a convective instability in which the 
source of free energy is the temperature gradient. In 
the more familiar Schwarschild convection, the entropy 
gradient is the source of free energy. For an arbitrar- 
ily oriented magnetic field, the disper sion relation of the 
MTI can be given bv IQuataertl (|2008h : 
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where TV is the Brunt- Vaisala frequency, 
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and 



Cj 2 = u? - (k ■ v A f , (4) 



where va — B/(47rp) 1 / 2 is the Alfven speed, and finally 
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(5) 



is the frequency for conduction to act on a given scale 
with b the unit vector directed along the magnetic field 
and x the thermal diffusivity 4 in units of cm 2 s~ 1 . This 
dispersion relation is written without loss of generality 

4 The literature is not consistent as to the use of \ an d k. This 
work will use \ to represent a true diffusion coefficient and n to 
represent a conductivity in units erg cm -1 s _1 K . 



for the geometry in which gravity and the initial atmo- 
spheric gradients are in the ^-direction and the initial 
magnetic field lies in the x-z plane. 

We now move to the limit of the MTI in which the 
temperature decreases with height in the atmosphere and 
the magnetic field is initially perpendicular to g. In this 
limit the growth rate of the instability simplifies to 



dlnT 

dz 



k 2 ' 



(6) 



which is growth at the isothermal limit of the Brunt- 
Vaisala frequency, where k± is defined as being per- 
pendicular to g. The MTI can be stabilized by strong 
magnetic fields, when magnetic tension begins to retard 
the convective motions, or by isotropic conduction that 
short-circu its the anisotropic driv ing n ecessary for this 
instabi lity. iParrish & Stond (|2005f ) and lParrish fc Stond 
(|2007al) describe local simulations of the growth and sat- 
uration of this instability in two- and three-dimensions. 
In short, the MTI saturates through a combination of 
three processees. First, the instability drives a magnetic 
dynamo that amplifies the field, increasing magnetic ten- 
sion. Second, if the system has adiabatic boundaries such 
that the temperature gradient can be decreased signifi- 
cantly, the system approaches isothermal and the driv- 
ing free energy term vanishes. Finally, in what perhaps 
is the most surprising result, the MTI reorients the di- 
rection of the magnetic field, preferentially aligning it 
with the background temperature gradients. This re- 
alignment reduces the growth rate of the instability as 
the field becomes more radial (aligned with the temper- 
ature gradient) and can be seen by examining equation 
(JSJ) ■ We shall see all of these saturation behaviors present 
in our study of the MTI in galaxy clusters. 

3. CLUSTER MODEL 

Before beginning any simulations of a cluster, the first 
step is to choose an appropriate model for the dark mat- 
ter and ICM that is inspired by observations and the- 
oretically consistent. Unfortunately, there is a bit of a 
"chicken and egg problem". What we are observing to- 
day is potentially already the saturated state that results 
from whatever instabilities and processes have been tak- 
ing place since the formation of the cluster. Nonetheless, 
we shall choose models that are generally consistent with 
the profiles observed today to study the general proper- 
ties of the ICM, as opposed to modeling any one specific 
cluster. 

Mu ch of th e inspi ration for our cluster model comes 
from iHughed (|1989h , in which the mass of the Coma 
Cluster is estimated from X-ray observations. Several 
ingredients go into the model which we shall discuss in 
turn: a profile for the dark matter, a profile for the gas 
properties, and a way to tie them together. Since major 
mergers are quite infrequent, the latter is quite easy to 
achieve, simply by assuming hydrostatic equilibrium. 

3.1. Dark Matter and Gas Profiles 

First, the dark matter profile: at the time of Hughes' 
paper, there was almost universal agreement that dark 
matter in clusters w as distribut ed according to the King 
dark matter profile (|Kind [l966). The profile is parame- 
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terized as 



Pdm(r) = PdmO 



1 



r 

Tdrr 



-i -n/2 



(7) 



where p,i m is the dark matter density binding the cluster, 
r dm is the scale radius, and n = 3, 4, or 5 in general. 
In the last decade simulations and some observational 
evidence have favored a somewhat cuspy central region 
given by the so-called NFW models (|Navarro et alll997f) 



Pdmjr) Sc 

Pent (r/r s )(l + r/r s y 



(8) 



where S c is the halo overdensity compared to the critical 
density at that epoch, and r s is the NFW scale radius. 
This work adopts the King model for the dark matter 
profile for several reasons. To begin with, we are not 
attempting to model the central region in detail, so the 
central region is less relevant. In addition, the cuspy 
dark matter profile makes it more difficult to hold the 
cluster in quiescent hydrostatic equilibrium. Finally, for 
the choice of n — 3 both profiles have pdm oc r -3 on 
scales larger than the core radius, and thus are roughly 
equivalent. This work will use the n = 3 King model for 
dark matter throughout. 

Again for the gas profile (the terms "gas" and 
"ICM" are used interchangeably in the literature) there 
are several possible choices. Hughes' paper makes 
the common choice of the iso thermal beta model 
(|Cavaliere fc Fusco-Femianolll976ft 



Pg( r ) = Pgo 



-3/3/2 



(9) 



where p g is the gas density, r c is t he gas scale rad ius, and 
the parameter is of order 0.75. iHughesI (|1989D derives 
cluster models by combining the King model with the 
isothermal beta model. There are two primary short- 
comings to this method. First, the resultant tempera- 
ture profiles are incredibly sensitive to the various free 
parameters, yielding many models that have unphysical 
temperatures that go to zero or infinity quite rapidly. 
Second, these models are in general not convectively sta- 
ble everywhere according to the Schwarzschild criterion, 
a poor initial condition for understanding a convective in- 
stability. It is also worth noting that these models arose 
during a time when X-ray measurements did not have 
sufficient resolution to resolve temperature gradients in 
clusters; thus, they are not necessarily as physically mo- 
tivated today. 

An appealing alternate model for the gas density is 
specified through a fit to the entropy profile. More 
recent X-ray observa tions are well fit by such profiles 
(|Donahue et al.ll2006f) . The general model is of the form 



S{r) = So 



(10) 



where S = pp 1 is the gas entropy, r c is the entropy scale 
radius, and a is a parameter. For a > 1, dS/dz > and 
the atmosphere is manifestly convectively stable. Out- 
side of the core many clusters are well fit observationally 
by a = 1.3, thus, that is the fiducial choice for this work. 
In general for the gas profiles there are three thermo- 
dynamic quantities: density, pressure, and temperature. 



One of these is determined by the equation of state, one 
of them is specified as an Ansatz, and the third must be 
calculated. When using the entropy model, the entropy 
is specified instead of one of these quantities. 

3.2. Generating Self- Consistent Cluster Models 

With appropriate choices of dark matter and gas pro- 
files, it is now time to tie them together through hydro- 
static equilibrium given as 



(11) 



where the gravitational potential, $, can be found by 
solving Poisson's equation, V 2< f> = AnGpdm, which in an 
isotropic spherical geometry is given as 



9$ _ 4ttG 
dr r 2 



dr'p dm (r')r n 



(12) 



Using the definition of the entropy as S = pp 7 , where 
7 is the adiabatic index, Eqn. |11| can be rewritten as 

dp = ldliig _ J_ p 2-^ 
dr 7 dr jS dr 

In order to simplify the subsequent derivations, we use 
the following dimensionless variables: p* = p/po, y = 
r/r dm , x = r/r c = yr dm /r c , where the subscript '0' 
refers to quantities at the cluster center. Solving Pois- 
son's equation for the gravitational acceleration, we ob- 
tain 

IT = 2 W)> ( 14 ) 

or 

where the function rj(y) is defined as 
rj{y) = arcsinh(y) — 



(15) 



Assuming the standard form for the pressure P = 
pksT J pmH where p, is the mean atomic number per nu- 
cleon (typically 0.6) and m# is the mass of the hydrogen 
atom, an equation for the density can be derived as 



dy 
dy 



a r drl 



C 



(2-7 



7 r c 1 + x a 



1 



y 



My), (16) 



where the coefficient C encodes many of the physical 
parameters of the cluster 



C 



_ AirGp dm opmHr 2 dri 



lk B T 

In more physical units C can be written as 



C = 15.0 



P_ 
0.6 



PdmO 



I'd,, 



2.0 x 10- 25 g cm- 3 

" '5/3\ /10keV\ 
~T) \ k B T ) 



(17) 



(18) 



500 kpc 

which is dimensionless. The equation for p can be for- 
mulated as either an ordinary differential equation or 
an integro-differential equation when the form of equa- 
tion Q2] is adopted for the gravitational acceleration. The 
equation must be solved numerically to calculate a pro- 
file; however, it is singular at y = 0. Equation 1161 can 
be integrated analytically in the limit of small y. This 
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TABLE 1 
Fiducial Cluster Parameters 



Parameter Name 


Symbol 


Value 


Dark matter scale radius 


^ dm 


750 kpc 


Entropy scale radius 


r c 


1000 kpc 


Central gas temperature 


To 


10 keV 


Dark Matter Mass 


M 


5.6 x 10 14 M 


Central dark matter density 


PdmO 


3.594 x 10" 26 g cm- 3 


Entropy profile exponent 


a 


1.3 


King exponent 


n 


3 


Constant for Eqn. 1161 


c 


6.064 



integration yields a non-singular starting point for inte- 
gration of 



P*( £ ) = l- 1 ( — 

1 \ r c 



(19) 



With this initial condition, the atmosphere is integrated 
using a fourth-o rder Runge-Kutta solver based on Nu- 
merical Recipes (|Press et al.lll986h . The fiducial param- 
eters for our cluster are given in Table [1] 

Figure [1] shows the calculated non-dimensional profiles 
for the fiducial cluster model detailed in this section. A 
nice property of this model is that the entropy is by con- 
struction increasing radially; thus, it is manifestly con- 
vectively stable by the Schwarschild criterion. 



500 1000 1500 
H(kpc) 



500 1000 1500 
H(kpc) 



500 1000 1500 

H(kpc) 



Fig. 1. — The non-dimensional temperature, density, and pres- 
sure profiles are plotted for the fiducial cluster model given by 
Table [T] 



4. METHOD 

4.1. Equations of MHD with Anisotropic Heat 
Conduction 

We solve the usual equations of ideal MHD with the 
addition of a vector heat flux, Q, and a gravitational 



acceleration g - 



djpv) 
dt 



dt 



V ■ (pv) = o, 





( S 2 > 




pvv + 







(20) 

pV$ = 0, 
(21) 



dE 
~dt 



D 2 



P- 



B(B-v) 



(22) 



VT / 4-7T 

+ V Q 

dB 

~dt 

where the symbols have their usual meaning. The total 
energy E is given as 

v • v B B 



/?V$ ■ v = 0, 
-Vx(»xB) = 0, 



(23) 



-3 



9 —1 

cm s . 



2 ■ 8* ' (24) 

and the internal energy, e = £>/(7 — 1). We assume 7 = 
5/3 throughout. 

The heat flux, Q = Q c + Q iso , where Q c , the 
anisotropic term is due to Coulombic conduction par- 
allel to magnetic field lines. The isotropic heat flux is 
given solely for comparison to the effect of anisotropic 
conductivity. The anisotropic conductivity is either set 
to a constant thermal diffusivity or given as the expres- 
sion 

Q c = -nk BX c(T, n)bb ■ VT, (25) 

where the diffusivity is given by the Spitzer conductivity 
(Spitzcr 1962) and b is a unit vector in the direction 
of the magnetic field. The Spitzer thermal diffusivity is 
parameterized appropriately for the ICM as 

Xc(T,n) = 8xl0 31 ( T \ (— 

,( 26 ) 

Note that this is a thermal diffusivity, and it is in- 
versely dependent on the density. When one speaks 
of the Spitzer conductivity, one often is speaking of 
K Sp — n ^BXc which has only the well-known T 5 / 2 de- 
pendence and no density dependence. 

4.2. Initial Conditions 

The initial conditions of the intracluster medium are 
described in detail in previous sections and the resulting 
plasma parameters are presented visually in Figure [1] A 
database from the Runge-Kutta solver is linearly interpo- 
lated to fill all space. The dark matter potential given by 
Equation [14] is applied as a static gravitational accelera- 
tion on the grid; thus, achieving hydrostatic equilibrium. 

We utilize two different magnetic field configurations 
in this study. The first initial magnetic field geometry is 
chosen to be both simple and compelling. Namely, the 
magnetic field is initialized to be purely azimuthal in the 
x-y plane at all radii such that r < r max . Essentially, 
this configuration is a sphere of nested circular magnetic 
field loops of constant strength. Numerically this config- 
uration is realized by utilizing the vector potential, where 
B = V x A. For a constant magnetic field of strength 
Bq defined as described, the vector potential is given as 

(27) 
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The vector potential is calculated at cell corners and then 
differenced to calculate the face-centered magnetic fields. 
The advantage of this magnetic field geometry is that the 
heat flux is zero in the initial configuration; thus, all sub- 
sequent evolution is due to the MTI. The fiducial mag- 
netic field strength is chosen to be 1 nG, possibly moder- 
ately higher than the poorly constrained primordial field 
but substantially less than today's observed field. 

A second magnetic field geometry is chosen to be a 
more realistic representation for the tangled fields ob- 
served in clusters today. There is little theoretical guid- 
ance for the magnetic field power spectrum in the ICM, 
so we choose a Kolmogorov power spectrum for the mag- 
netic field, initialized in Fourier space as 



A(k) = A 



vpeak 



(28) 



where fc poa k is chosen as the wavenumber corresponding 
to 2-4 times the grid scale. We utilize the Fast Fourier 
Transform (FFT) as was done in lRoettiger et alj (jl999l ) 
to calculate the vector potentials in real space. We also 
randomize the phase to avoid correlating the modes. The 
magnetic energy scales as B 2 oc fc -21 " -1 ). We initialize 
each Gaussian component separately, so for the compo- 
nentwise Kolmogorov spectrum, the appropriate choice 
is a — —17/6 to give the familiar scaling for energy of 
fc-U/3_ T hi s fc-space scaling for Kolmogorov turbulence 
is appropriate for initializing the ID power spectrum (as 
opposed to the 3D power spectrum) in each direction 
separately. We numerically check that the divergence 
of the magnetic field remains zero and renormalize the 
magnetic energy such that (B 2 ) = (B 2 ), 

4.3. Timescales 

In this section we discuss several of the charac- 
teristic timescal e s for a cluster as was outlined in 
iParrish fc Stond (|2007al ). The Brunt- Vaisala frequency 
can be written in a more useful form for the cluster prob- 
lem as 

N2 _ i^m»s 

7 or or 

The maximum growth rate is given by the isothermal 
limit of the Brunt- Vaisala frequency, namely 

Table [2] examines these timescales in the galaxy cluster 
at a fiducial radius of 450 kpc. Several of these timescales 
are dependent on the length scale of interest. Thus, they 
are considered on global scales, roughly a wavelength of 
1 Mpc, and the scale of the mean free path, roughly 30 
kpc. The magnetic field is assumed to be 1 nG for these 
calculations. 

As can be seen by examination, the MTI has tens of 
growth times during a Hubble time, allowing significant 
rearrangement of the atmosphere. The magnetic field 
estimated using primordial values clearly plays no role, 
even on the length scale of the mean free path. 

On mean free path timescales, the conduction time is 
faster than the sound crossing time. The potential worry 
is due to th e thought that the heat flux is saturated in 
the sense of ICowie fe McKed (|1977l ). The right compar- 
ison there is not the conduction time to the mode sound 



crossing time but the electron thermal velocity to the ve- 
locity needed for electrons to transport the heat. For the 
atmosphere considered here, the saturated heat flux is 
several orders of magnitude larger than the actual heat 
flux. 

4.4. Numerical Methods 

We use the 3D vers i on of the Athena M HD code 
(jGardiner fc Stond[2008t [St one fc Gardinerl2008h to sim- 
ulate the MTI far into the nonlinear regime. Both the 
MHD and the heat transport methodolo gy and tests have 
been d escribed in these references and IParrish fc Stone! 
(2005), respectively The thermal diffusivity is either set 
as a constant or permitted to vary spatially according 
to the standard Spi tzer prescription. Here we a dopt the 
adopt the method of Sharma 'fc Hammettl (|2007h for han- 
dling variable thermal conductivities. Namely, the prod- 
uct of density and thermal diffusivity is interpolated to 
the cell faces using a harmonic average, e.g. in one di- 



mension as 



("X)»+i/2 (nx)i (nx)i 



(31) 



This harmonic averaging prevents the Courant condition 
from becoming severe due to discontinuous diffusivities 
and densities at interfaces, should they develop. The 
timcstcp is determined with respect to the maximum 
thermal diffusivity on the grid and the conduction mod- 
ule is sub-cycled. 

All the runs described in this chapter are performed on 
a uniform Cartesian three-dimensional grid with the ini- 
tial cluster atmosphere described in Sj3] The domain ex- 
tends from -3.0 x 10 24 cm to +3.0 x 10 24 cm (968 kpc) in 
each direction. Most of the runs presented in this paper 
are simulated at a resolution of (128) 3 or corresponding 
to a zone size of 15 kpc respectively. The simulations are 
typically run for 3 x 10 17 s or 9.5 Gyr, a substantial frac- 
tion of cosmic time. The magnetic field is initialized in a 
region r max which is slightly less than the length of the 
cube. Typically for our simulations r max = 2.5 x 10 24 
(806 kpc). Thus, the initial magnetic field is confined 
to a spherical region. Since there is no magnetic field 
outside of r max , there is initially no conduction in this 
region. 

We use modified reflecting boundary conditions for all 
the MHD variables, in which the pressure and density 
are extrapolated in the ghost zones. In order to prevent 
any negative values for pressure or density in the ghost 
zones, we introduce pressure and density floors applied 
only in the ghost zones. The magnetic field and veloc- 
ity field components are simply reflected at the bound- 
ary. The heat fluxes are forced to be adiabatic by setting 
Q ■ N = on the boundary. 

To seed multiple modes of the MTI and break sym- 
metry, we add Gaussian white noise perturbations to the 
velocity field such that the applied perturbation is a fixed 
fraction of the local sound speed everywhere within the 
cluster. A typical perturbation is 1% of the sound sound 
speed in local velocity magnitude. In the absence of 
applied perturbations, it is reassuring to note excellent 
Cartesian symmetry is maintained. 

4.5. ICM Atmosphere Stability for MHD-only Runs 
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TABLE 2 

TlMESCALES OF THE ICM 



Timescale 


Definition 


A = 1 Mpc (Myrs) 


A = 30 kpc (Myrs) 


Sound crossing time 


T 3 


A 


610 


18 


Alfven crossing time 


TA 


_ c i 


5.4 X 10 7 


1.6 x 10 6 


Thermal conduction time 


T x 


— x 


2.3 X 10 3 


2.1 






xc. 






Brunt-Vaisala time 


TJV : 


= N- 1 


810 


810 


Growth time 


Ttr - 


-l 

- Gmax 


880 


880 



In order to both show that these numerical methods are 
valid and to provide a basis for comparison, we introduce 
a reference calculation in which the atmosphere is initial- 
ized as previously described, except thermal conduction 
is disabled. In this case, the MTI is not active and the 
stability properties of the background atmosphere can be 
deduced. 

The background state is run for for a time of 3 x 10 17 
s (9.6 Gyr) on a (64) 3 grid, at which point it is analyzed. 
The cluster remains relatively quiescent. Of course, there 
is a slight mismatch resulting from discretizing a spher- 
ical problem on a Cartesian grid. The maximum Mach 
number reached anywhere on the grid is 9.8 x 10 -3 and 
the rms Mach number is even smaller. Another charac- 
terization of the quiescence is the change in the central 
pressure. Over the entire run time the maximum change 
in the central pressure from the initial condition is only 
0.2%. It is also useful to look at the conservation of 
energy. It is not just internal energy that must be con- 
sidered, but also the total gravitational binding energy 
of the cluster: 



~t 



1 

e H — pv ■ v 



1 



-B ■ B + p$ 



0. 



(32) 



where the angle brackets reflect a volume average. For 
the case where the MTI is active, the density profile is 
modified, and the last term is important to consider. By 
this measure, the total energy is conserved to better than 
0.1% for both the MHD-only and MHD plus conduction 
cases. 

Finally, it is useful to introduce a diagnostic at this 
point which reflects the net alignment of the magnetic 
field with the radial direction. The diagnostic is mathe- 
matically represented as (\b ■ r|), which is the overlap 
of the magnetic field unit vector and the radial unit 
vector, averaged over the cluster region in which the 
field was originally initialized. The original azimuthal 
field construction would have this quantity be ideally 
zero, but discretization errors give an initial value of 
(\b ■ r|) = 2.3 x 10~ 4 for (64) 3 . At the conclusion of 
the MHD only run, this value has risen to 0.02, a small 
change to be compared to the runs with conduction later. 
In addition, in the MHD-only run, there is approximately 
a 5% decay in magnetic field energy through reconnec- 
tion, which is also the primary source of the change of 
the magnetic field orientation. Thus, in this section, 
we have characterized the evolution of the cluster at- 
mosphere with MHD-only runs, and shown that this at- 
mosphere remains very quiescent. 

5. SIMULATIONS OF THE MTI 



In order to simulate the effects of the MTI on the intr- 
acluster medium, We have performed a number of sim- 
ulations that are listed in Table [31 The key runs of in- 
terest are the ones with purely anisotropic conductivity 
labeled A1-A8. Run Al has a fixed thermal diffusivity 
and provides the cleanest physics laboratory for under- 
standing the MTI. Runs A2-A7 have a variable Spitzer 
thermal diffusivity. These runs also constitute a resolu- 
tion and geometry study. Run A8 is the same as run A2, 
except the magnetic field is three orders of magnitude 
larger. For variable conductivity, the global conduction 
timestep is limited by the largest conductivity anywhere 
on the grid, which has made (256) 3 achievable only with 
several tens of thousands of processor hours. 

We also provide a number of runs for comparison pur- 
poses. First, the fiducial run we have already outlined in 
§4.51 is labeled F0. This run shows what the background 
state is in the absence of any thermal conduction. We 
also provide runs 11—12 which utilize purely isotropic con- 
ductivity at a fixed fraction of the Spitzer conductivity 
to highlight the difference between the instability and 
simple isotropic conductions. 

5.1. Fixed Anisotropic Thermal Conduction 

We begin our discussion with the simplest illustrative 
case of this paper, Run Al. This simulation of the 
MTI was performed with a fixed thermal diffusivity of 
8.0 x 10 31 cm 2 s _1 . In addition, the magnetic geome- 
try is initialized with simple circular magnetic field lines. 
Figure [2] shows the dramatic effect of the MTI on the ra- 
dial shell-averaged temperature profiles. In a time short 
compared to the Hubble time, the temperature profile is 
quite dramatically rearranged. By 7.9 Gyr, the initial 
temperature profile is almost completely erased, and the 
cluster driven towards isothermality. Concomitantly the 
density in the central regions is driven somewhat higher 
with the central density peaking at 6.6 x 10 -27 gcm~ 3 . 
While initially some heat is transported inwards, Brehm- 
strahlung cooling scales like n 2 ; thus, with the addition 
of cooling, this density increase may enhance the cooling 
flow problem. 

Simultaneously, the MTI is driving a magnetic dynamo 
which acts to amplify the magnetic field. This is quanti- 
fied as 

where the angle brackets denote a volume average over 
the sphere with r < r max . For Run Al, the magnetic 
energy density is amplified by a factor of 18.4 from the 
initial energy density. There is actually somewhat more 
magnetic dynamo action, as some of the amplified flux is 
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TABLE 3 

Initial Properties of Nonlinear Runs 



Run 


Conduction Type 


Conductivity (cm 2 s _1 ) Field Geometry B / (4tt) 1/2 (fiG) 


Resolution 


FO 


None 





Circular 


lO" 9 


64 3 


AI 


Anisotropic 


8 x 10 31 


Circular 


lO" 9 


128 3 


A2 


Anisotropic 


Spitzcr 


Circular 


lO" 9 


128 3 


A3 


Anisotropic 


Spitzer 


Tangled 


10 -9 


128 3 


A4 


Anisotropic 


Spitzcr 


Circular 


lO" 9 


64 3 


A5 


Anisotropic 


Spitzer 


Circular 


lO" 9 


256 3 


A6 


Anisotropic 


Spitzer 


Tangled 


lO" 9 


64 3 


A7 


Anisotropic 


Spitzer 


Tangled 


lO" 9 


256 3 


A8 


Anisotropic 


Spitzer 


Circular 


10" 6 


128 3 


11 


Isotropic 


8.0 X 10 30 T 5 / 2 n -1 


Circular 


lO" 9 


128 3 


12 


Isotropic 


2.4 X 10 31 T 5 / 2 n" 1 


Circular 


10" 9 


128 3 






TABLE 4 












Saturated Properties of Nonlinear Runs 








Run a/r2\ 
rtun o \r> ) 


RMS Mach a Max. Mach a 


(b ■ f) /s P itz b 


*max.KE (Gyr) 






FO 0.95 


2.4 x 10 -3 9.8 x 10~ 3 


0.05 









AI 18.4 


0.03 0.3 


0.62 0.28 


7.7 






A2 19.0 


0.04 0.3 


0.59 0.28 


7.0 






A3 17.2 


0.03 0.1 


0.67 0.45 


1.8 






A4 13.2 


0.04 0.3 


0.60 0.43 


5.2 






A5 23.5 


0.04 0.3 


0.59 0.26 


7.1 






A6 14.6 


0.05 0.2 


0.70 0.49 


1.7 






A7 59.8 


0.05 0.3 


0.68 0.43 


1.8 






A8 0.95 


0.04 0.2 


0.44 0.13 


7.7 






11 1.1 


1.0 x 10~ 3 1.2 x 10~ 3 


0.19 


1.0 






12 1.1 


1.1 x 10~ 3 1.4 x 10~ 3 


0.19 


0.9 





a Evaluated at the peak of kinetic energy. 

b Evaluated at the peak of ([&• r|), normalized to the instantaneous Spitzer heat flux. 



saturation properties of this run and all other simula- 
tions are given in Table [U This magnetic field increase is 
significant, but clearly not entirely sufficient to raise the 
primordial field to its present-day value. Certainly, much 
of the increase comes from flux-freezing during gravita- 
tional collapse, which is not included here. In addition, 
smaller-scale field increases could come from purely ki- 
netic processes such a s the Alfven wave cascade or ot her 
kinetic process (e.g., iSchekochihin fc Cowlevl l2?K)6affbh . 
The large scale turbulence can serve as energy input for 
this process. 

5.2. Spitzer Anisotropic Conduction 

We proceed to discuss the more realistic fiducial case 
of run A2. Here, instead of a fixed conductivity as in AI, 
we utilize the full variable Spitzer conductivity as given 
in Eqn. |26| . We begin by examining the thermal evolu- 
tion of A2 in Figure [31 The evolution of the temperature 
profile is somewhat smoother than in the case with fixed 
conductivity and has approached a fairly isothermal clus- 
ter by about 7 Gyr after the start of the run. Just as a 
reminder, the anisotropic conduction is initially zero due 
to B = beyond r max , which is slightly greater than 800 
kpc. 

Figure |4] shows the time evolution of the volume- 
averaged magnetic and kinetic energy for run A2. The 
magnetic energy is continually amplified during the en- 
tire evolution, reaching a peak amplification of 5{B 2 ) rs 
19.0. The kinetic energy reaches a peak around 7.7 
Gyr. At this time, the RMS Mach number, where 
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Fig. 2. — The MTI dramatically rearranges the radial tempera- 
ture profiles for Run AI. 



transported by penetrative convection to r > r max and 
is not included in this average amplification factor. The 
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Fig. 3. — The shell-averaged radial temperature for circular 
magnetic fields with Spitzer conductivity in run A2 evolves signif- 
icantly. 

Mach = H/c s , averaged over the cluster volume is 0.04. 
The maximum Mach number of the turbulence at any 
given zone in the cluster is 0.3. The kinetic energy de- 
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Fig. 4. — The volume- averaged evolution of the magnetic and 
kinetic energy during the course of run A2 show magnetic field 
amplification and significant motions. 

cays somewhat as the source of free energy, the temper- 
ature gradient, has been mostly exhausted. In practice, 
clusters are constantly being heated by gravitational in- 
fall and AGN activity in the cluster core. These heat- 



ing terms are neglected but clearly would contribute to 
maintaining the temperature gradient and source of free 
energy; thus, the decay of turbulence is likely overpre- 
dicted in the simulation. 

Perhaps the most interesting part of the evolution dy- 
namics is that the magnetic field geometry is significantly 
rearranged by the MTI. Initially, the magnetic field ge- 
ometry was almost perfectly azimuthal, lying completely 
in the x-y plane. At the conclusion of the simulation, 
the magnetic field is highly biased in the radial direc- 
tion, with (|6 • r|) = 0.62 at maximum. Figure [5] shows 
the evolution of this metric. The tendency of the MTI 
to align the magnetic field with the background tempera- 
ture gradient was previo usly seen in the 3D stratifi ed box 
simulations presented in iParrish &: Stonel (|2007bf ). This 
large scale rearrangment of the magnetic field has poten- 
tial observational consequences which will be discussed 
later. 



0.6 - 
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Fig. 5. — The volume-averaged radial component of the magnetic 
field, (\b- r\) , is amplified significantly during the course of run A2. 

It is notable to compare this simulation to the fiducial 
simulation, FO, in which conduction does not play a role. 
In the absence of a dynamo driven by conduction, the 
magnetic field actually decays by reconnection. In FO, 
the flow remains very quiescent and the magnetic field 
almost completely retains its initial azimuthal geometry. 

The same simulation as run A2 has also been per- 
formed at (64) 3 (run A4) to provide a resolution study. 
In general, as is to be expected, the prime difference is 
the degree of amplification of the magnetic field. Small 
scale wavelengths grow fastest for the MTI; thus, higher 
resolution permits more small-scale field growth at some- 
what faster timescales. The higher resolution also per- 
mits greater tangling of the field by convective motions. 
The Mach number is slightly higher in lower resolution 
runs; whereas the magnetic field geometry remains fairly 
constant. Given the dyanmically weak magnetic fields 
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achieved after 9.5 Gyr of growth, we are not yet in the 
limit where magnetic tension is playing a significant role 
on the small scales. Thus, the magnetic field amplifi- 
cation is not yet saturated with the tangled field am- 
plification resulting in S(B 2 ) « 59.8 at 256 3 and 19.0 
at 128 3 . Higher resolution, while not computational 
tractable, will likely yield even higher magnetic field am- 
plification. Aside from this exception, we are not seeing 
any significant trends with resolution. 

5.3. Heat Flux Diagnostics 

We also analyze the cluster heat flux from our stored 
time slices at the conclusion of a run. We will focus on the 
radial component of the heat fluxes in a spherical shell of 
width 100 kpc, centered on a radius of 450 kpc from the 
cluster center. We begin by defining a fiducial conduc- 
tivity to be the Spitzer conductivity in this shell if the 
conduction were isotropic at the Spitzer value, namely: 

dT 



Q r = -nk B xc(T,n) 



dr 



(34) 



This value is the same as for anisotropic thermal con- 
duction with purely radial field lines. In Figure [6l we 
compare the heat flux due to conduction to this Spitzer 
value. However, there remains some ambiguity as the 
temperature profile evolves over the course of the simu- 
lation. Thus, we normalize the conductive heat flux to 
both the initial Spitzer heat flux and the instantaneous 
Spitzer heat flux. We define the Spitzer fraction, the 
oft-quoted / paramater as 

/Spitz = Qcond/Q- (35) 

In Table IU we utilize the instantaneous fiducial heat flux 
at the time of maximum (|6 • f|)for normalization of this 
quantity. Figure [6] shows that as the ICM temperature 
gradient is depressed, the conductive heat flux normal- 
ized to the instantaneous heat flux increases substan- 
tially, until /spitz is larger than 0.4. 

Meanwhile, calculating the heat flux due to the trans- 
port of enthalpy is worth a careful treatment. The ad- 
vection of enthalpy can be written as 



Q 



enthalpy 



7 



7-1 



(36) 



where we neglect the transport of kinetic and magnetic 
energy since the flows are both low mach number and 
high-/3, respectively. We can split this heat flux into an 
advective term proportional to a mass flux (ex nv) and 
a term that is a purely convective heat flux. Writing the 
pressure in Eqn.[3B] as nk B T, we define several quanti- 
ties: 



n = (n) + 5n 
v = (v) + Sv 
T=(T)+6T, 



(37) 
(38) 
(39) 



where, for example (n) is the spatial average of the num- 
ber density in the shell and 6n is the local deviation of 
a quantity from its average. In this formulation the new 
total heat flux becomes 

7 



Q r = 7 - ± ^k B (nTv r ) (40) 
-k B (((n) + 5n)((T) + ST)({v r ) + Sv r )' 
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Fig. 6. — Here we show the evolution of the radial heat fluxes 
for run A2 in a 100 kpc shell located 450 kpc from the cluster 
center. The conductive heat flux normalized to the initial fiducial 
heat flux (black, solid line) reaches a maximum before decreasing, 
while normalized to the instantaneous fiducial heat flux (red, long- 
dashed line) continues to increase. The heat fluxes due to mass 
advection (blue, dotted line) and convection (green, short-dashed 
line) are shown normalized to the initial fiducial heat flux. 

In future equations, we drop the r subscript. In simpli- 
fying this expression, terms that are proportional to the 
average of a delta quantity, like (8n) vanish by definition. 
We can then separate it into two quantities: one due to 
a mass flux, identified as terms that scale as nv, and the 
remaining terms into a convective flux. Thus, 

7 -k B ((n)(T)(v) + (T)(5nSv)), 



adv ' 



Qc 



7-1 

7 



k B {(v)(8nST) + (n)(6v6T) 



(41) 
(42) 



7-1 
+ (SnSTSv)) , 

where the second term of equation (f42|) dominates the 
convective flux. 

Figure [6] shows these two heat fluxes normalized to 
the initial fiducial heat flux. The mass advection heat 
flux is enhanced due to the inability of the dark matter 
profile to adjust to changes in the ICM structure, so it is 
not a physically important quantity. What is interesting, 
however, is that the convective heat flux is quite small. In 
fact, it is roughly four orders of magnitude smaller than 
the conductive heat flux. The small convective heat flux 
seems to result from two principal causes. First, unlike 
in solar convection, there is a second channel open for 
energy flow, namely directly through conduction which 
is quite efficient. Second, as the magnetic fields become 
more radial, the buoyant driving is significantly reduced, 
reducing the baseline convective motions. 

At this point it is also appropriate to discuss run A8. 
Run A8 is chosen to be initialized at the MTI stability 
boundary for strong magnetic fields, namely kvA ~ c max , 
for wavenumbers corresponding to scales slightly smaller 
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than the cluster radius. In this limit, all but the largest 
scale modes are suppressed entirely by magnetic tension. 
There is both theoretical motivation and numerical ev- 
idence that the conductive flux increases as the initial 
magnetic field increases. We do indeed find that for run 
A8, the convective flux is within an order of magnitude of 
the conductive flux; however, this phase lasts only a very 
short time (less than 500 Myr) and the net convective flux 
never is larger than 0.035 Q. In addition, for this run, the 
total heat flux carried is considerably smaller than the 
other runs. Thus, we conclude that conduction is the 
dominant heat transport mechanism and that convective 
heat transport can never be larger than a small fraction 
of the Spitzer conductivity for MTI-driven turbulence in 
galaxy clusters. 

5.4. Comparison to Isotropic Conductivity 

In order to emphasize the nature of the MTI as 
a truly convective instability, we have performed sim- 
ulations with purely isotropic conductivity. A long- 
favored approach to understanding thermal conduction 
in galaxy clusters has been to assume a tangled mag- 
netic field geometry passively determined by hydrody- 
namic turbulence. This tangled geometry then pro- 
vides an effective therm al conductivity that is a spec- 
ified fraction of Spitzer (iRechester &: Rosenbluth] 119781 : 
IChandran fe Cowlevlll998t iNaravan fc Medvedevli200lD 
A major point of this paper is that such an approach is 
incorrect since thermal conduction along the magnetic 
field self-consistently modifies the magnetic field geome- 
try, resulting in an /s p itzor that evolves with time and is 
not consistent with these theories. Thus, here we present 
simulations II and 12 which have a thermal diffusivity 
fixed at 10% or 30% of the Spitzer value, respectively, 
for comparison purposes. 

The temperature profiles of these two runs evolve in 
ways quite similar to that of run A2 and A3, reaching 
fairly isothermal in 5-7 Gyr. A much more illustrative 
difference is found by considering the time evolution of 
the two runs. The runs are initialized with a Mach 0.01 
velocity perturbation that decays in time quite quickly, 
peaking at 1 Gyr or less, as opposed to 7 Gyr for run A2. 
Without the MTI, the magnetic field is barely amplified 
at all, S(B 2 ) w 1.3 for both cases. 

Finally, the magnetic field geometry is somewhat rear- 
ranged by the cluster motions, but does not even reach 
a geometrically isotropic state, with (\b ■ r\) reaching a 
maximum of 0.19. Of course, we started in a state of 
pure azimuthal fields, so the key message here is that a 
field that started isotropically distributed would remain 
isotropically distributed — there is no force here pushing 
the field to be radially biased as with the MTI. Clearly, 
the comparison between these isotropic conduction cal- 
culations and the preceeding anisotropic heat conduction 
calculation shows that the MTI plays a significant role in 
setting a self-consistent solution for magnetic field ampli- 
fication, convective velocities, and for the magnetic field 
geometry of a cluster. 

5.5. Tangled Magnetic Fields 

We now turn to the case of tangled magnetic fields, ex- 
amined in run A3. This run has the exact same set-up as 
the previously discussed runs, except the intial magnetic 



field is a turbulent magnetic field distributed according 
to the Kolmogorov scaling law. Thus, in this initial case, 
there is a zeroth order anisotropic heat flux that exists 
independent of the MTI's evolution. Figure [7] shows the 
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Fig. 7. — The evolution of the shell-averaged radial temperature 
for tangled magnetic fields with Spitzer conductivity in run A3 is 
quite rapid. Note that the times plotted are roughly a factor of 
two earlier than for run A2 in Figure [3] 

fairly rapid evolution of the temperature profile for the 
tangled field, reflecting this zeroth order radial heat flux. 
The cluster reaches a fairly isothermal temperature in 
3.5 Gyr. While the temperature gradient is being sup- 
pressed, the magnetic energy is amplified by a modest 
factor of 17.2, excluding the flux pumped to radii larger 
than r max . Figure [8] shows the evolution of the magnetic 
and kinetic energy over the run. The kinetic energy peaks 
at 1.8 Gyr, much earlier and slightly higher in amplitude 
than the circular field case. The short wavelength-modes 
that exist in the tangled field geometry grow and saturate 
on somewhat faster timescales than the largest modes in 
the problem. 

An interesting property of the resultant turbulence is 
seen in examining Figure [U The kinetic energy at its 
peak is significantly larger than the magnetic energy. At 
the peak of the kinetic energy, the kinetic energy has been 
amplified by a factor of 72.9. The Alfven Mach number, 
Ma = (v/va) can be simply determined as the square 
root of the ratio of the kinetic to magnetic energy. At 
the peak of the kinetic energy, the turbulence is highly 
super- Alfvenic, with Ma ~ 230.. This property appears 
to be generally true as long as magnetic tension is not 
dominant. 

It is quite intriguing that the re-orientation of the mag- 
netic field geometry can be observed quite vividly even 
for the tangled field case. Starting with an initial geome- 
try with (|6 • r|)= 0.5, the radial component is amplified 
to a peak value of 0.67 at 2.1 Gyr, before decaying after 
the free energy is exhausted. Figure [9] shows the evolu- 
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Fig. 8. — The volume- averaged evolution of the magnetic and 
kinetic energy during the course of run A3. 

tion of this quantity over the simulation. Simultaneously, 
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Fig. 9. — The volume-averaged radial component of the magnetic 
field, {\b ■ r|), is amplified even for the tangled magnetic field in 
run A3. 

we examine the evolution of the heat fluxes in this prob- 
lem in the same way we did for run A2. Clearly, the 
initial conductive heat flux is already a modest fraction 
of Spitzer due to the isotropically tangled magnetic field 
lines. As the simulation evolves, the conductive heat flux 
normalized to the instaneous fiducial heat flux reaches a 
value of /spitz = 0.45 at the peak of • f\). The radial 



component continues to rise further past 0.5 as the tem- 
perature gradient is further depleted in this simulation. 
Again the convective heat fluxes are quite small in both 



or 
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Fig. 10. — This figure demonstrates the evolution of the radial 
heat fluxes for run A3 in a 100 kpc shell located 450 kpc from the 
cluster center. The conductive heat flux normalized to the initial 
fiducial heat flux (black, solid line) reaches a maximum before 
decreasing, while normalized to the instantaneous fiducial heat flux 
(red, long-dashed line) continues to increase. The heat fluxes due 
to mass advection (blue, dotted line) and convection (green, short- 
dashed line) are shown normalized to the initial fiducial heat flux. 

an instantaneous and time-average sense. These simula- 
tions in a tangled field geometry prove the robustness of 
the MTI in tangled fields when the power is initially on 
a variety of scales. 

6. OBSERVATIONAL CONSEQUENCES AND 
DISCUSSION 

We have shown in the previous section that properly 
accounting for anisotropic thermal conduction in the in- 
tracluster medium drives unique physics with potential 
observational consequences. The resultant physics is sig- 
nificantly different from the pure MHD case without con- 
duction or the isotropic conduction case. Of course, with- 
out heating or cooling mechanisms included, the simula- 
tions presented here are insufficient to make predictions 
for any particular cluster. Instead, we briefly outline here 
two potential observational consequences that will be a 
focus of future work. 

6.1. Turbulent Velocities 

In this work, we predict rms turbulent velocities of 
order 3-4% the speed of sound and peak velocities of 30- 
40% of the sound speed. In a relaxed cluster that has not 
undergone a major merger, this may be the dominant ve- 
locity component in the ICM; although, certainly wakes 
from passing galaxies may also be a significant contribu- 
tion. The typical ICM velocities for the MTI are pre- 
dicted to be approximately an order of magnitude larger 
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than the velocities that would result from pure isotropic 
conduction alone. 

To determine the measurability of these hydrodynamic 
motions, we must compare the Doppler broadening due 
to hydrodynamic motions to the simple thermal Doppler 
broadening. The most success will come in measuring 
hydrogen- and helium-like iron lines since the thermal 
broadening is much smaller for the heavy iron nuclei. 
We can express the ratio of the broadening mechanisms 



Ai/ t 



urbulcnt 



Afthermal 



= 0.46 



(Vrms \ 
60 kms- 1 / 



k B T 
5keV 



-1/2 



1/2 



56 vtih 
( 43 ) 

These two contributions are roughly equal for an rms 
turbulent velocity of 60 km s _1 and a temperature of 
1 keV. For the peak velocity fluctuations, the hydrody- 
namic motions dominate; however, at the rms level, the 
velocity fluctuations will be harder to see. Clever model- 
ing of the actual line shapes, as in llnogamov &; Sunvaevi 
(2003), may improve the detectability threshold of hy- 
drodynamic turbulence. 

As of yet, there have been no successful measurements 
of ICM velocities. Had the X-ray calorimeter on Astro- 
E2 (Suzaku) not failed, it would have been able to image 
the velocity fluctuations of the ICM at the level of the 
thermal or hydrodynamic broadening. The next chance 
for confirming any of these predictions will occur with 
Constellation- X. Its spectral resolving power is energy- 
dependent, but based on current performance require- 
ments, should be able to resolve the predicted level of 
velocity fluctuations in iron lines with energies around 6 
keV. 

6.2. Magnetic Field Geometry 

Several authors have noted a discrepancy in the mea- 
surements of cluster magnetic fields, finding the fields de- 
rived from rotation measure (RM) observations are typi- 
cally a factor of four to ten higher than fi elds derived from 
inverse Compton (IC ) measurements ()Petrosia n1 [200l 
iCarilli &~T avlor 20Q2J). The efforts at explaining this dis- 
crepancy are not particularly satisfactory. 

First, consider fields derived from a rotation measure. 
The RM can be written as 



RM = 812 



» B ■ dl radians m 



(44) 



where n e is given in c m 3 , B is measured in //G, and 
dl is measured in kpc (ICarilli fc Tavlorll2002ft . The ro- 
tation measure is proportional to the component of the 
magnetic field along the line of sight times the electron 
density. Typically, the field geometry is assumed to be 
isotropic in direction. An enhancement in the radial field 
along the line of sight would represent an enhancement 
in rotation measure, and an overestimate of the magnetic 
field strength. 

Alternatively, the magnetic field may be derived from 
inverse Compton emission from the CMB off relativistic 
electrons. While the inverse Compton emission has no di- 
rect dependence on the magnetic field, the magnetic field 
strength may nonetheless be derived by comparing the 
inverse Compton emission to the synchrotron emission, 
provided they are both produced by the same population 
of relativistic electrons. More illumination is provided by 



examining the formulae for the luminosity of both phc 
nomena: 

4 2 2 

L sync = -p 7 co-tUb, 



Lic = ^P 7 c(T T U 1 , 



(45) 
(46) 



where (3 — v/c, 7 = (1 — /3 2 ) -1 / 2 , <tt is the Thompson 
scattering cross section, Ub = B 2 /8tt is the magnetic 
energy density, and f7~ is the energy density of the pho- 
ton field. The similarity of these two equations is quite 
evident, giving 



Lie 



J sync 



U 1 
Ub~' 



(47) 



The inverse Compton emission peaks in the X-ray around 
20 keV, corresponding to electrons with a 7 ~ 5000. The 
corresponding sy nchrotron emission peaks in the radio 
around 100 MHz (jBagchi et al.lll998h . 

The IC emission is entirely independent of the mag- 
netic field geometry; however, the synchrotron emission 
is a weak function of magnetic field geometry. Clearly, 
more detailed analysis will be required in the future; how- 
ever, the basic idea is intriguing. By simply positing a 
magnetic field that is preferentially radial to the line of 
sight, the rotation measure-derived field is overestimated 
while the IC/synchrotron-derived magnetic field is not 
significantly changed. Future work will examine whether 
the action of the MTI can explain this observational dis- 
crepancy. 

6.3. Temperature Profiles 

It may seem as if the temperature profile is the most 
obvious observerable, however we have a "chicken and 
egg" problem. Our initial conditions posited in JO] re- 
flect our observations of clusters today. Thus what we 
are observing are relaxed clusters that could be in the 
saturated state of any instability that has taken place. 
In fact, the initial conditions for a relaxed cluster today 
should probably be taken from a typical disturbed clus- 
ter (e.g. Perseus) that has recently undergone a major 
merger. Additionally, we have at this point neglected 
heating sources due to minor mergers, AGN activity, 
and infall that can be important in heating a cluster. 
Nonetheless, we find that significant evolution of the tem- 
perature profile can occur in the relevant time between 
cluster major merger events. In future work, we will 
work to more systematically connect our predicted tem- 
perature profiles with observations. 

6.4. Thermal Evaporation and Missing Baryons 

There have been several recent suggestions that 
Sunyaev-Zel'dovich (SZ) and X-ray measurements of 
galaxy clusters are incompatible. As mentioned before, 
the X-ray measurements are sensitive to the line inte- 
gral of the square of the density; whereas, the SZ signal 
is directly related to the line integral of t he pressure (or 
therm al energy) of the baryons. Recently. Afshordi et al.l 
(2007) have used composite WMAP data to find that 
35 ± 8% of the cluster baryons are missing from the 
ICM. Namely, the authors argue that the gas fraction 
of clusters is consistently poor relative to the universe 
as a whole. The validity of this measurement will likely 
be confirmed or refuted in the near future as much more 
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accurate SZ surveys, including the Atacama Cosmology 
Telescope (ACT) and the South Pole Telescope (SPT) as 
well as numerous others, begin to make maps and mea- 
surements of clusters. 

This piece of observational evidence has motivated pro- 
posals for th ermal e vapor ation of galaxy clusters, such as 
the work by lLoebl (|2007f h This work supposes that the 
outer portions of galaxy clusters could thermally evapo- 
rate, losing up to 10% of their mass. The primary pro- 
posed mechanism for this is leakage of suprathermal par- 
ticles with large Coulomb mean free paths. While this 
work is perhaps not theoretically robust, it does require 
magnetic fields that are preferentially radial in direction 
in order for this large mean free path argument to trans- 
late into suprathermal particles escaping. Radial mag- 
netic fields are of course one of the predictions of the non- 
linear MTI simulations to date, thus satisfying one key 
requirement for this process. Reconsidering this prob- 
lem with a methodology similar to Parker's treatment of 
the solar wind could place this process on a more firm 
theoretical footing. Of course, observations and theory 
are not yet well-developed on this topic, but it is worth 
keeping these facts in mind until the observations become 
more certain. 

7. SUMMARY AND FUTURE WORK 

The intracluster medium is a dilute magnetized plasma 
in which the mean free path is a moderate fraction of 
the size of a cluster. In this plasma, heat flows exclu- 
sively along magnetic field lines making the ICM suscep- 
tible to the magnetothermal instability in regions where 
the temperature profile is decreasing outwards. We have 
performed fully three-dimensional time-dependent MHD 
simulations of the intracluster medium and outlined its 
evolution and potential observational consequences. We 
have also performed simulations without conduction and 
with only isotropic conduction to provide a basis for com- 
parison and understanding of the MTI. 

There are several common features of the evolution of 
the MTI in the intracluster medium. First, the MTI 
is capable of rearranging the temperature profile of the 
ICM on timescales significantly shorter than the Hub- 
ble time. Second, the MTI drives subsonic convection 
that amplifies the magnetic field in a dynamo. Third, 
the magnetic field geometry is driven to a highly-radially 
biased field. Finally, we have self-consistently calculated 
/spitzer ^ 0.4 — 0.5 for the first time and find that conduc- 
tion can be very efficient at transporting energy. Mean- 
while, convective heat transport by the MTI in clusters is 
not efficient at all. These conclusions simultaneously help 
understand several puzzles surrounding the ICM includ- 
ing the general flatness of the large-scale temperature 
profiles and the origin of the present-day cluster mag- 
netic fields. 

We have also discussed several potential observational 
consequences of the MTI and anisotropic transport. 
First, future X-ray satellites such as Constellation- X 
may be able to image the subsonic velocity fluctuations. 
Second, future work may explain the discrepancy be- 
tween inverse Compton and rotation measure determina- 
tions of the magnetic field strength. Finally, a better un- 
derstanding of the masses of galaxy clusters may be ob- 
tained by understanding the ICM physics in much more 
detail. Careful measurements of galaxy cluster masses 



can provide detailed cosmological tests. 

An important lesson from this work is that high-/3 
magnetic fields, though weak in terms of pressure, are not 
at all negligible. In the long mean free path regime, such 
as in the intracluster medium, the physics of anisotropic 
transport along magnetic field lines can self-consistently 
drive the temperature profile evolution, magnetic field 
amplification and magnetic field geometry, as well as tur- 
bulence. Treating the magnetic field as a simple passive 
contaminant in hydrodynamic turbulence misses key rel- 
evant physics. 

The future work here is manifold. In the very near 
future, we will shift our attention inward to the cluster 
cores. This region of the cluster, in which the temper- 
ature increases outward, is unstable to a cousin of the 
MTI kno wn as the heat-flux-driven buoyancy insta bility, 
or HBI (jQuataertl 120081 : iParrish fc QuataertJl2008l ). In 
contrast to the MTI, this instability drives field lines to 
be perpendicular to the background temperature gradi- 
ent, potentially exacerbating the cooling flow problem. 
In the short term, we will study the effects of the HBI 
and cooling in cluster cores. For the global MTI cal- 
culations in particular, several steps remain. First, a 
study of the cores of clusters including heating and cool- 
ing processes will likely require static mesh refinement. 
Second, we have yet to do a thorough parameter study. 
In such a study, we would quantify the variation of the 
MTI physics as parameters such as the cluster mass and 
central temperature are varied. This type of understand- 
ing is very important for the use of clusters as standard 
candles in cosmology. Finally, we are very interested in 
exploring the observational consequences of the MTI in 
more detail. In particular, we would like to understand 
theoretically the modification of the magnetic field mea- 
surement through inverse Compton and rotation measure 
methods, as the geometry of the magnetic field changes 
from isotropic to radially-biased. With a theoretical un- 
derstanding, we could compute exact IC and rotation 
measure maps from a high-resolution simulated cluster 
for comparison with observational data. 

In short, clusters of galaxies provide a rich physics lab- 
oratory for testing Braginskii-MHD, the MTI, and many 
other phenomena. We have developed analytical and 
computational tools for understanding these phenomena, 
and we will continue to apply these extended MHD mod- 
els to improving our understanding of clusters of galaxies. 
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